A degenerate three-level laser with a parametric amplifier 



Eyob Alebache"wEI and K. Fesseha 
Department of Physics, Addis Ababa University, P. O. Box 33085, Addis Ababa, Ethiopia 

(Dated: February 9, 2008) 

The aim of this paper is to study the squeezing and statistical properties of the light produced 
by a degenerate three-level laser whose cavity contains a degenerate parametric amplifier. In this 
quantum optical system the top and bottom levels of the three-level atoms injected into the laser 
cavity are coupled by the pump mode emerging from the parametric amplifier. For a linear gain 
coefficient of 100 and for a cavity damping constant of 0.8, the maximum intracavity squeezing is 
found at steady state and at threshold to be 93%. 
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I. INTRODUCTION 

There has been a considerable interest in the analysis 
of the quantum properties of the squeezed light gener- 
ated by various quantum optical systems [l], [3,11, 0, 
H, IS B Hi • In squeezed light the fluctuations in one 
quadrature is below the vacuum level at the expense of 
enhanced fluctuations in the other quadrature, with the 
product of the uncertainties in the two quadratures satis- 
fying the uncertainty relation. In addition to exhibiting 
a nonclassical feature, squeezed light has potential appli- 
cations in precision measurements and noiseless commu- 
nications ncum. 

Some authors have studied the squeezing and statisti- 
cal properties of the light produced by three-level lasers 
when either the atoms are initially prepared in a coherent 
superposition of the top and bottom levels [HI, [IH, [l4| 
or when these levels are coupled by a strong coherent 
light [13j. These studies show that a three- level laser 
can under certain conditions generate squeezed light. In 
such a laser, three-level atoms in a cascade configuration 
are injected at a constant rate into the cavity coupled 
to a vacuum reservoir via a single-port mirror. When 
a three-level atom makes a transition from the top to 
bottom level via the intermediate level, two photons are 
generated. The two photons are highly correlated and 
this correlation is responsible for the squeezing of the 
light produced by a three-level laser. On the other hand, 
it is well known that a parametric oscillator is a typical 
source of squeezed light [3, [3|, |J, [5j, [fj , with a maximum 
intracavity squeezing of 50%. Recently Fesseha [13] has 
studied a three-level laser with a parametric amplifier in 
which three-level atoms, initially prepared in a coherent 
superposition of the top and bottom levels, are injected 
into the cavity. He has found that the effect of the para- 
metric amplifier is to increase the intracavity squeezing 
by a maximum of 50%. 

In this paper we consider a degenerate three- level laser 
whose cavity contains a degenerate parametric amplifier 



(DPA) and coupled to a vacuum reservoir. The top and 
bottom levels of the three-level atoms injected into the 
cavity are coupled by the pump mode emerging from the 
parametric amplifier. And the three-level atoms are ini- 
tially prepared in such a way that the probabilities of 
finding the atoms at the top and bottom levels are equal. 
We expect that a highly squeezed light can be gener- 
ated by the quantum optical system under consideration. 
Thus our interest is to analyze the squeezing and statis- 
tical properties of the light generated by this system. 

We obtain, applying the master equation, stochastic 
differential equations for the cavity mode variables asso- 
ciated with the normal ordering. The solutions of the 
resulting equations are used to determine the quadrature 
variance, the squeezing spectrum, and the mean photon 
number. Moreover, applying the same solutions, we de- 
termine the antinormally ordered characteristic function 
with the aid of which the Q function is obtained. Then 
the Q function is used to calculate the photon number 
distribution. 



II. STOCHASTIC DIFFERENTIAL EQUATIONS 

Three-level atoms in a cascade configuration are in- 
jected into the laser cavity at a constant rate r a and 
removed from the cavity after a certain time r. We rep- 
resent the top, middle, and bottom levels of a three-level 
atom by |a), |6) and |c), respectively. We assume the 
transitions between levels \a) and \b) and between levels 
| b) and |c) to be dipole allowed, with direct transitions 
between levels \a) and |c) to be dipole forbidden. We 
consider the case for which the cavity mode is at reso- 
nance with the two transitions * \b) and \b)— > |c) (see 
Fig. 1). 

The Hamiltonian describing the coupling of levels \a) 
and |c) by the pump mode emerging from the parametric 
amplifier can be expressed as 



H' = i^(\c){a\-\a){c\), 



(1) 



'Electronic address: yob a@yahoo.com 



in which f2 = 2<?'/i with g' and /i being respectively the 
coupling constant and the amplitude of the pump mode. 
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FIG. 1: A degenerate three-level laser with a degenerate para- 
metric amplifier. 



In addition, the interaction of a three-level atom with the 
cavity mode can be described by the Hamiltonian 

H" = ig[tf(\b){a\ + \c)(b\) - 6(|a)(6| + |&)(c|)], (2) 

where g is the coupling constant and a is the annihila- 
tion operator for the cavity mode. Thus the Hamiltonian 
describing the interaction of a three-level atom with the 
cavity mode and with the pump mode emerging from the 
parametric amplifier has the form 



ff = i 5 [at(|6)( a | + |c>(6|)-o(| a )(6| + |6)(c|; 
+ ij(\c)(a\-\a)(c\). 



(3) 



We take the initial state of a single three-level atom to 



be 



1^(0)) 



1 



1 



V2 V2 

and hence the density operator for a single atom is 



(4) 



p A (Q) = ~\a)(a\ + \\a){c\ + \\c){a\ + l -\c){c\. (5) 

It can be readily established that the equation of evo- 
lution of the density operator for the laser cavity mode, 
coupled to a vacuum reservoir, has in the linear and adi- 
abatic approximation the form [l5j 



dt 



p = R{2d) pa — da) p — pad)) 

+ S{2dpd) — a) hp — pa) a) 

+ U{d) par + apa — pa) 2 — a 2 p) 



+ V(d) pa) + apa — pa 1 — a) J p), 



t2, 



where 



R = 



A 
4~B 



(6) 



(7a) 



U 



V 



A 
4B 

A 

is 
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2 2 
2 2 



2 
2 



B = (l+f3 2 )(l + ^), 

f3 = n/ 7 , 



A 



r 



(7c) 

(7d) 

(7e) 
(7f) 

(8) 



is the linear gain coefficient, k is the cavity damping con- 
stant, and 7 is the atomic decay rate assumed to be the 
same for all the three levels. 

Moreover, a degenerate parametric amplifier with the 
pump mode treated classically is describable in the inter- 
action picture by the Hamiltonian 



IE 



# = ~(at a -a'), 



(9) 



in which e = Xp with A being the coupling constant. The 
master equation associated with this Hamiltonian has the 
form 



d » e i 
dt P= 2 if,a 



d 2 p + a) 2 p - pa) 2 ) 



ta-i 



(10) 



Now on account of Eqs. (J6J) and (|T0|) . the master equation 
for the cavity mode of the quantum optical system under 
consideration can be written as 



-j t P = 2 ~a p + a) p- pa) ) 
+ R(2d) pa — ad) p — paa)) 
+ S(2apa) — a) dp — pa) a) 
+ U (a) pa) + apa — pa) 2 — a 2 p) 
+ V(a) pa) + apa — pa 2 — a) 2 p). 



(11) 



We next proceed to determine, applying this master 
equation, the stochastic differential equations for the cav- 
ity mode variables. To this end, applying (fTTj) one readily 
finds 

±(a) = (R-S){a) + (U-V + e){tf), (12) 



dt 



(a 2 ) = 2{R-S){a 2 )+2{U-V + e){d)a)+e-2V, (13) 



S ~ AB 



2nB 
~A~ 



(7b) 



dt 



(a T a) = 2{R-S)(a^d) + {U-V + s){(d u ) + (a 2 )) + 2R. 

(14) 



We note that these equations are in the normal order and 
the corresponding c-number equations are 



dt 



(a) = -(S-R)(a) + (U-V + e)(a*), (15) 



— (a 2 ) = -2(5 - R)(a 2 ) + 2(U - V + e)(a*a) + e - 2V, 
dt 

(16) 



(a* a) = -2(S-R)(a*a) 

+ (U - V + e){(a* 2 ) + (a 2 )) + 2R. (17) 



dt 



On the basis of Eq. (fTS]) , one can write the stochastic 
differential equation 

ja(t) = -(S-R)a(t) + (U-V + e)a*(t) + f(t), (18) 

where f(t) is a noise force the properties of which remain 
to be determined. We observe that Eq. (fTS")) and the 
expectation value of Eq. (fT8|) will have the same form if 

(/(*)) = 0. (19) 
Moreover, it can be readily verified using (fT8| that 
d 



dt 



(a 2 (t))=-2(S-R)(a 2 (t)} 

+ 2(U -V + e)(a*(t)a{t)) 
+ 2(a(t)f(t)), (20) 



and 



-{a*(t)a{t))=-2(S-R){a*(t)a(t)) 

+ (U-V + s)((a* 2 (t)) + (a 2 (t))) 
+ (a(t)f*(t)) + (a*(t)f(t)). (21) 

Comparison of Eqs. (fTS"]) and (|20"jl as well as Eqs. (fT7|) 
and l|2ip shows that 

{a{t)f{t)) = \{e-2V), (22) 



(a(t)f*(t)) + (a*(t)f(t)) = 2R. (23) 

Furthermore, a formal solution of Eq. (fT5|) can be writ- 
ten as 

a(t) = tt (0) e -( s - fi " 

+ [ e -(s-R)(t-t') [([/ _ F + e )a*(0 + /(*')]*'• 
Jo 

(24) 

Using Eq. along with (I22p . one easily finds 



e-< s - fl X t - t '></(f)/(f))df = 5(e - W). (25) 
o z 



Based on this result, one can write [12,15] 

(f(t')f(t)) = (e-2V)S(t-t'). (26) 



It can also be established in a similar manner that 



</(*)/*(*')> =2R6(t-t'). 



(27) 



We note that Eqs. (|26|) and (|27|) describe the correla- 
tion properties of the noise force f(t) associated with the 
normal ordering. 

Now introducing a new variable defined by 

a±(t) = a*(t)±a(t), (28) 

we easily get with the help of (18) that 

d 



dt 



a±(t) = -\ T a±(t) + f*(t)±f(t), (29) 



where 

\ T = (S-R)T(U-V + e). (30) 
The solution of Eq. (|29|) can be written as 



a±W=a±(0)e- A ^ + fe-^( M '»(r(t')±/(0)*' 
Jo 

(31) 

It then follows that 

a(t) = A(t)a(0) + B(t)a*(0) + F(t), (32a) 



in which 



and 



A(t) = i( e - A -' + e- A +'), 



B(t) = ~(e- x - t -e- x +% 



F(t)=F+(t)+F_(t), 



(32b) 



(32c) 



(33a) 



with 



F±(t) = i jT e-W-'Hrtf) ± /(tO)*'. (33b) 



III. QUADRATURE FLUCTUATIONS 

In this section we seek to calculate the quadrature vari- 
ance and squeezing spectrum for the cavity mode under 
consideration. 
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A. Quadrature variance 



The variance of the quadrature operators 
a + = a) + a 

and 

a_ = i(a' — a) 



(34) 



(35) 



is expressible in terms of c-number variables associated 
with the normal ordering as 



A<4 



l±{a±(t),a±(t)), 



(36) 



in which a± (t) is given by Eq. (|28p . Assuming the cavity 
mode to be initially in a vacuum state and taking into 
account (j3Tj) together with (fl"9|) . wee see that 



(«±(t))=0. 
In view of this result, Eq. ([36]) reduces to 

A4 = i±(o4(t)). 

Furthermore, employing Eq. (f29|) . one easily gets 

|<4(*)> = -2\ T {al(t))+2{a ± (t)r(t)) 
±2(a ± (t)f(t)). 



(37) 
(38) 



(39) 



With the aid of Eq. along with (J2BJ) and (J27K, we 

readily obtain 

j f (a 2 ± (t)) = -2A T (a 2 ± (t))+2{e-2V± 2R) . (40) 

The steady-state solution of this equation turns out to 
be 



(a\(t)) 



(41) 



e-2V±2R 



Now on account of gTj) together with fl30j), Eq. lj38| 
takes at steady state the form 



Aa 



= 2k(1 + /3 2 )(1 + /3 2 /4) + A(A + /3 2 ) 

~ (2k - 4e)(l + /? 2 )(1 + /3 2 /4) + A(2(3 - p) 

(42a) 



and 



Aa' 



2 K (l + /3 2 )(l + /3 2 /4) + 3A/3 2 



(2k + 4e)(l + /? 2 )(1 + /3 2 /4) + A(4/3 + /3 3 ) ' 

(42b) 

where we have used Eqs. ([7a|) - ([7e|) . 

In the absence of the parametric amplifier the quan- 
tum optical system under consideration reduces to just 
a degenerate three-level laser. The quadrature variance 
for this system has upon setting (3 = e = in (|42a[) and 
(I42bp the form 



2A 



(43a) 




FIG. 2: Plots of the quadrature variance [Eq. ()46b|l ] versus 
j3 for k = 0.8 and A=100 in the presence of the pump mode 
and in the absence of the nonlinear crystal (dotted curve) 
and [Eq. ()51b|l ] for k = 0.8 and A=100 in the presence of the 
parametric amplifier (solid curve). 




FIG. 3: Plots of the quadrature variance [Eq. (|51b|) ] versus (3 
for n — 0.8 and for different values of the linear gain coefficient 
A in the presence of the parametric amplifier. 



and 



Aai = 1. 



(43b) 



Since neither of the quadrature variance is less than one, 
the light produced by the degenerate three-level laser 
is not in a squeezed state. We therefore observe that 
the particular initial preparation of the three-level atoms 
we have considered does not lead to the generation of 
squeezed light. We next consider the case in which three- 
level atoms are not injected into the cavity. Thus setting 
A = (3 = (with (i ^ 0) in Eqs. ([I2"a|) and (|42b) . we 
have 



Aa 2 



2s 



(44a) 



5 



and 



Aa 2 



2e 



(44b) 



At threshold, e — k/2, the above equations reduce to 

Aa+ — > oo (45a) 



and 



(45b) 



The squeezing in this case is exclusively due to the para- 
metric amplifier. 

On the other hand, when the nonlinear crystal of the 
parametric amplifier is removed from the cavity, the 
quantum optical system under consideration reduces to 
a coherently driven three-level laser. Hence upon setting 
e = Xfi = (with fi ^ 0) in Eqs. (|4"2a|l and (142bj) the 
quadrature variance for this system becomes 



Aa 2 



and 



Aa 2 = 



2 K (l + /3 2 )(l + / 3 2 /4) + A(4 + /3 2 ) 
2k(1 + (3 2 ){l + /3 2 /4) + A{2{3 - 03) 



2k(1 + /3 2 )(1 + /3 2 /4) + 3A/3 2 
2k(1 + (3 2 ){l + /3 2 /4) + A(4/3 + /3 3 ) ' 



(46a) 



(46b) 



It is interesting to consider the case for which the coher- 
ent light driving the three-level laser is sufficiently strong. 
Thus we note that for /3 = O/7 3> 1, Eq. (|46ap reduces 
to 



Aa% 



f + A/p 2 



(47) 



Since 1 is very large, we can drop the term A/ f3 2 in Eq. 
(1471) , so that 



Aai 



K-2A//3 

Following the same procedure, we also get 
Aa 2 _ = 



■2A/0 



(48) 



(49) 



Now comparison of Eqs. (|44a[) and (|4"8")) as well as Eqs. 
(|44b|) and (|49|) shows that a degenerate three-level laser 
driven by a strong coherent light behaves in exactly the 
same manner as a degenerate parametric oscillator, which 
is in agreement with the assertion made in Ref. [l6l ]. 
However, as can be seen clearly from Fig. 2, a relatively 
high degree of squeezing occurs for small values of (3. 

Inspection of Eq. ([30)) shows that A+ is nonnegative 
while A_ can be positive, negative or zero. We thus note 
that Eq. (|29[) will not have a well-behaved solution if 
A_ < 0. Now setting A_ = and taking into account 
([30j). we get 



e = — 



A(2f3-f3 3 ) 



4(1 + /? 2 )(1 + /3 2 /4)' 



(50) 



We can then interpret this as the threshold condition for 
the system under consideration. Therefore expressions 
(|42ajl and (|42bj) take at threshold the form 



and 



Aa 2 _ 



2k(1 + /? 2 )(1+/3 2 /4) + 3A/3 2 
4k(1 + /3 2 )(1+/3 2 /4) + 6A/3 



(51a) 



(51b) 



The dotted curve in Fig. 2 represents the quadrature 
variance for the coherently driven three-level laser and 
the solid curve in the same figure represents the quadra- 
ture variance for the degenerate three-level laser with the 
parametric amplifier at threshold. It is easy to see from 
Fig. 2 that the presence of the parametric amplifier leads 
to better squeezing. Moreover, the minimum value of the 
quadrature variance, described by Eq. l(5Tb|l . for A=100 
and k = 0.8 is found to be Aa 2 _ = 0.068 and occurs at 
(3 = 0.067. This result implies that the maximum intra- 
cavity squeezing for the above values of the linear gain 
coefficient and cavity damping constant is 93% below the 
vacuum level. In addition Fig. 3 shows that the degree 
of squeezing increases with the linear gain coefficient. 



B. Squeezing spectrum 

The squeezing spectrum of a single-mode light is ex- 
pressible in terms of c-number variables associated with 
the normal ordering as 

/■oo 

S° ± ut (oj) = l±2Re (a ^(t), a ^(t + r)) ss e^ T dT, 
Jo 

(52a) 

where 



t (t) = a* out (t)±a out (t). 



(52b) 



For a cavity mode coupled to a vacuum reservoir, the 
output and intracavity variables are related by 



c4 u *(t) = v^o±(t). 



(53) 



Now taking in to account Eqs. (|37|) and (|53|> . the squeez- 
ing spectrum can be put in the form 

/•OO 

Sg*(u) = l±2KRe / (a±(t)a±(t + T)) ss e luJT dT. (54) 
Jo 



On the other hand, the solution of Eq. ([29)1 can also be 
written as 



a±{t + r) = a±(t)e- x * T + [ , 

Jo 

±f{t + T'))dr'. 



-A T (r-r') 



(/*(* + t') 



(55) 



Upon multiplying ([55)) by a± (t) and taking the expecta- 
tion value of the resulting equation, one readily obtains 
at steady state 



(a±(t)a±(t + t)) ss = (a±(t)) ss e~ 



(56) 
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FIG. 4: Plots of the squeezing spectrum [Eq. (|57b|l ] versus /3 
for ft = 0.8 and A=25 in the presence of the pump mode and 
in the absence of the nonlinear crystal (dotted curve) and [Eq. 
(|58b|) ] for k = 0.8 and A=25 in the presence of the parametric 
amplifier (solid curve). 
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FIG. 5: Plots of the mean photon number [Eq. (|61[) ] at steady 
state versus f3 for k = 0.8, and A=25 in the presence of the 
pump mode and in the absence of the nonlinear crystal (dot- 
ted curve) and for k = 0.8, and A=25 in the presence of the 
parametric amplifier with e = 0.3 (solid curve). 



Now using Eqs. (|4T|) and ([56)) , the squeezing spectrum is 
found to be 



ST'M = 1 + 



2n{e-2V + 2R) 
X 2 _ + uj 2 



(57a) 



and 



S „„, M = 1 _M£_1^). m 

A + + UJ" 

Applying Eqs. ([SO)). (ffa |) -(|7e ]l . and ([30)). the squeezing 
spectrum ([57)) can be put at threshold in the form 



S° + Ut {u) = 1 + 



, A(4+/3 2 ) 
~ 2 



and 



S°_ ut (uj) = 1 



,4/3(3-3/3/2) ^ 

(i +/ a2)(i+/3 2 /4)' 



r )2 +W 2- 



(58a) 



(58b) 



(l+/3 2 )(2+/3V2)' 



Inspection of this equation shows that there is perfect 
squeezing for /3 = uj = and for any values of A and k. 
The dotted curve in Fig. 4 shows that there is almost 
perfect squeezing at /3 = 0.016 and the solid curve indi- 
cates that the presence of the parametric amplifier loads 
to perfect squeezing at j3 = 0. 



IV. PHOTON STATISTICS 



A. The mean photon number 



Applying Eq. ([32ap and its complex conjugate, the 
mean photon number of the cavity mode, assumed to 
be initially in a vacuum state, can be written as 



= (F*(t)F(t)). 



(59) 



On account of Eqs. (|33ap and l|33b[) together with the 
correlation properties of the noise force f(t), we readily 
obtain 



(F*(t)F(t)) = 2R -|K±£ (i_ e -aA-«) 



4A_ 
2R + 2V -e 



(l-e- 2X+t ). (60) 



4Aj 



In view of this result, the mean photon number takes the 
form 

, 2e(l + (3 2 )(l + l3 2 /4) + A(2-[3 + (3 2 /2 + /3 3 /2) 
^ a; 4(1 + (3 2 )(1 + (3 2 /4)(k - 2e) + 2A(2/3 - /3 3 ) 
x (l-e- 2A -*) 

2e(l + /3 2 )(1 + /3 2 /4) + ^(2/3 - 3/3 2 /2 + /3 3 /2) 
4(1 + /3 2 )(1 + /3 2 /4)(k + 2e) + 2,4(4/3 + /3 3 ) 
x(l-e^ 2A+t ). (61) 



We now proceed to calculate the mean photon number 
and the photon number distribution for the cavity under 
consideration. 



Fig. 5 clearly indicates that the parametric amplifier 
contributes significantly to the mean photon number for 
relatively small values of (3. 
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B. The Photon number distribution 

We finally seek to calculate, employing the Q function, 
the photon number distribution for the cavity mode. The 
photon number distribution for a single-mode light is ex- 
pressible in terms of the Q function as [TtJ 



P(n,t) 



2n 



n\ da* n da 1 



-[Q(a*,a,t)e Q * a ] 



(62) 



Now using (f62l) and (jAllj) . the photon number distribu- 
tion for the cavity mode can be written in the form 



P(n,t) 



2n 



_ (c 2 - d 2 ) 1 ' 2 d 

nl da* n da n 
x exp[(l - c)a*a - d(a* 2 + a 2 )/2] a , =a=0 - (63) 



Upon expanding the exponential functions in power se- 
ries, we have 



P(n,t) 



(c 2 - d 2 ) 1 ' 2 ^ (-l)'+ m (l - c) k d l+r 



n 

Q2n 



E 

kirn 



2 l + m k\l\m\ 



da* n da r 



{a* k+2l a k+2l } a * =a=Q , 



(64) 



so that on carrying out the differentiations and applying 
the condition a* = a = 0, there follows 



P(n,t) 



{c 2 -d 2 y/ 2 (-i) i+m (i - c ) k d l+m {k + 2iy. 



k! in 



2 l + m k\l\m\{k + 21 - n)\ 



(k + 2m)\ 
(k + 2m — n)\ 



(65) 



Applying the properties of the Kronecker delta symbol 
and the fact that a factorial is defined for nonnegative 
integers, we obtain 



P(n,t) = (c 2 -d 2 )^ 2 ynl 



(1 



-21 j2l 



2=0 



' 2 2l l\ 2 {n - 2l)V 



(66) 



where [n] = n/2 for even n and [n] = (n — l)/2 for odd 
n. As can be seen from Fig. 6, the steady-state photon 
number distribution decreases with the photon number. 
Moreover, the probability of finding even number of pho- 
tons is in general greater than the probability of finding 
odd number of photons. Although the photons are gen- 
erated in pairs in this quantum optical system, there is 
a finite probability to find odd number of photons in- 
side the cavity. This is because some photons leave the 
cavity through the port mirror. It also appears that the 
presence of the parametric amplifier increases the prob- 
ability of finding even number of photons and decreases 
the probability of finding odd number of photons. 



V. CONCLUSION 

In this paper we have considered a degenerate three- 
level laser whose cavity contains a degenerate parametric 




FIG. 6: Plots of the photon number distribution [Eq. (|66[l ] at 
steady state versus photon number for A=100, /3 — 0.067, and 
k = 0.8 in the presence of the parametric amplifier with e = 
0.3 (solid curve), and e = in the presence of the pump mode 
and in the absence of the nonlinear crystal (dotted curve). 



amplifier, with the top and bottom levels of the three- 
level atoms coupled by the pump mode emerging from 
the parametric amplifier. We have obtained using the 
master equation stochastic differential equations. Ap- 
plying the solutions of the resulting equations, we have 
calculated the quadrature variance and squeezing spec- 
trum. Moreover, using the same solutions, we have de- 
termined the mean photon number and the photon num- 
ber distribution. We have found that the light generated 
by the quantum optical system is in a squeezed state, 
with the maximum intracavity squeezing attainable be- 
ing 93%. The parametric amplifier increases the squeez- 
ing significantly over and above the squeezing achievable 
due to the coupling of the top and bottom levels by the 
pump mode. In addition, we have seen that there is per- 
fect squeezing of the output light for j3 = uj = and 
for any values of A and n. We have also found that the 
presence of the parametric amplifier leads to a significant 
increase in the mean photon number for small values of /?. 
The plots of the photon number distribution show that 
the probability of finding odd number of photons is in 
general less than the probability of finding even number 
of photons. Moreover, the same plots indicate that the 
probability of finding even number of photons is greater 
in the presence of the parametric amplifier than in the 
absence of the nonlinear crystal. However, the opposite 
of this assertion holds for the probability of finding odd 
number of photons. 
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APPENDIX A: THE Q FUNCTION 



We evaluate the Q function applying the relation 

Q(a*, a, t) = [ dz 2 (j>(z*,z,t)exp(z*a — za*), (Al) 

7T ./ 



where 



(A2) 



is the antinormally ordered characteristic function de- 
fined in the Hiesenberg picture. Using the identity 



(A3) 



this function can be written in terms of c-number vari- 
ables associated with the normal ordering as 

<f>{z*, z, t) = e- z * z (exp(za*(t) - z*a{t))). (A4) 

Since (|15[) is a linear differential equation, we see that 
a(t) is a Gaussian variable [HI]. In addition, on account 
of (|32a[) . (|33a|) . (|33b[) . and (JT9J) we easily find the mean 
value of a(t) to be zero. Hence a(t) is a Gaussian variable 
with vanishing mean. One can then express (|A4j) in the 
form [l|| 

4>(z*,Z,t) = e~ z ' z exp(^([z 2 a* 2 + z* 2 a 2 - 2zz*aa*})). 

(A5) 

Now employing Eqs. (|32a[) . (|33a|) . and (|33b|) along with 
the correlation properties of the noise force f(t), one can 
establish that 



4A_ 



(A6) 



(a* a) = 
+ 



2fl-2F + £ 2A 

4A_ 1 J 
2R + 2V~e 



4A, 



(1 _ e -2A +t) _ (A7) 



Therefore, with the aid of Eqs. (|A6j) and (|A7|), Eq. (|A5 
can be put in the form 



(z*,z, t) = exp(-az*z - b{z z + z* 2 )/2), (A8) 



in which 



, 2R-2V + e^ 2A t . 
a = l + — — (1-e-^-*) 

2R + 2V-e^ _ 2A+tN 



4A, 



(A9) 



b= 2E-2^ + £(1 _ e _ 2A _ t) 



4A_ 
2R + 2V -e 
4A~ 



(l- e - 2A +*) 



(A10) 



Upon substituting (|A8| into (|A1|) and performing the 
integration, the Q function is found to be 

f c 2 _ (f\l/2 d 

Q(a* , a, t) — exp[— ca*a — —(a* 2 + a 2 )), 

7T 2 

(AH) 

where 

(A12) 



a 2 - b 2 ' 



d = 



a 2 -b 2 ' 



(A13) 
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